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Abstract 

The reduced Bethe-Salpeter equation with scalar confinement and 
vector gluon exchange is applied to quark-antiquark bound states. 
The so called intrinsic flaw of Salpeter equation with static scalar 
confinement is investigated. The notorious problem of narrow level 
spacings is found to be remedied by taking into consideration the 
retardation effect of scalar confinement. Good fit for the mass spec- 
trum of both heavy and light quarkomium states is then obtained. 

PACS number(s): 



I Introduction 



To understand quark confinement is the most important task in QCD and hadron 
physics. Lattice QCD calculations show that the interquark potential for a heavy quark- 
antiquark pair QQ in the static limit is well described by a linear confining potential, plus 
a short-ranged Coulomb potential [HJ. Phenomenologically, these potentials have been used 
with the Schrodinger equation for nonrelativistic heavy quarkonium systems like cc and 
bb states, and satisfactory results for their mass spectrum have been obtained. Not only 
the spin-independent but also the spin-dependent QQ potentials are studied both in lattice 
QCD|| and the quark potential model[|3], f|, ||. Most results seem to be consistent with the 
picture that the dominant part of linear confinement potential is transformed as a Lorentz 
scalar, while the Coulomb potential stems from one gluon exchange, which has the feature 
of a Lorentz vector. In particular, the fact that the spin-orbit term (Thomas procession 
term) induced by the scalar confining potential tends to partially compensate the spin-orbit 
term generated by one gluon exchange, has been strongly supported by the observed fine 
splittings of P-wave cc and bb states p, [|, ||. 

However, a rather serious problem seems to remain if the spin-independent relativistic 
correction, caused by the static (instantaneous) scalar confining potential, is taken into 
consideration. This spin-independent term is of the same order as the Thomas precession 
term in the nonrelativistic expansion in terms of p 2 jvn 2 . As noted before @, including this 
spin-independent term 
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Hsi = - jW* + *Sf + -f- + ^|), (1) 
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where S(r) is the static scalar confining potential and is usually assumed to take the form 

of S(r) = Ar with A being the string tension, into the Hamiltonian will badly disturb the 

mass spectrum of mesons (even for cc states), because this term is negative and unreasonably 

large for higher excited states, making the level spacings for higher lying states unreasonably 

small. This problem is probably due to the fact that the scalar confining potential has been 

treated as an instantaneous potential, which is valid in the static limit but may not hold 

when relativistic corrections are taken into account. This problem was also noted by other 

authors 1§ in the framework of the reduced Salpeter equation. This equation is equivalent 

to the Breit equation to the first order of p^/m 2 , and may be used to study higher order 

relativistic corrections for systems containing the charm quark and even lighter quarks. It 

was found || that in the framework of the reduced Salpeter equation with an instantaneous 

scalar confining potential, the level spacings (e.g., the 25 — IS spacing) would tend to vanish 

for qq mesons when the constituent quark mass approachs to zero, and difficulties are already 



evident for the cc states. It was then pointed out|| that there is an intrinsic flaw in the 
approach which uses the reduced Salpeter equation with static scalar confinement potential. 

To overcome this difficulty, several scenarios have been put forward |||7]]. The chief 
differences between those works are of the useage of interaction potentials. By now, there are 
no mature theories or calculations for qq confinement interaction in QCD, and customarily 
used potentials are phenomenological and have some uncertain perameters in them. For 
different procedure of evaluating these parameters, one can have some different ways to 
fit the experimental data. Therefore, to classify some of the most effective alternatives of 
quark- ant iquark interaction potentials seems still premature. 

Despite of the limited understanding for confinement at present, more theoretical efforts 
should be made to study this problem. In our opinion, the difficulty with the reduced Salpeter 
equation and the static scalar confinement is probably due to the improper treatment that 
the confining interactiton is purely instantaneous. With some retardation effect of quark 
confinement being considered, even within the framework of reduced Salpeter equation the 
level spacings for qq mesons could become normal since the retardation effect might cancel 
the sick disturbance caused by the spin- independent correction from the instantaneous part 
of scalar confinement || . To implement this idea, we will assume the confinement kernel in 
momentum space to take the form 

G(g)oc 1 1 
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where q is the 4-momentum exchanged between the quark and antiquark in a meson. In fact, 
this form was suggested for the dressed gluon propagator at small momenta to implement 
confinement 0] . Here we will use the same form but regard it as an effective scalar confine- 
ment kernel. Then, if the system is not highly relativistic we may make the approximation 

and may further express go i n terms of its on-shell values which are obtained by assuming 
that quarks are on their mass shells. This should be a good approximation for cc and bb 
states, because they are nonrelativistic systems and the binding energies are smaller than 
the quark masses, therefore the quarks are nearly on their mass shells. In order to get a 
qualitative feeling about the retardation effect considered here, we will also use ([3]) for light 
quark mesons, though the approximations are not as good as for heavy quark mesons. With 
above approximations, the scalar confinement kernel becomes instantaneous again but we 
have incorperated some retardation effect into the kernel. In the static limit, the retardation 
term vanishes and the kernel returns to G(q) oc rap, which is just the Fourier transform of 
the linear confining potential. 



In this paper, we will use this modified scalar confining potential in which the retardation 
effect is incorperated and the one-gluon-exchange potential in the framework of the reduced 
Salpeter equation to study the mass spectra of qq mesons including both heavy and light 
mesons. We will concentrate on the 0~ and 1~ mesons to examine their level spacings. 

II Reduced Salpeter equation with scalar and vector 
interactions 

In quantum field theory, a basic description for the bound states is the Bethe-Salpeter 



equation ||10 1 . Define the Bethe-Salpeter wave function of the bonud stase | P) of a quark 



ip(xi) and an antiquark ip(x 2 ) as 

X (x 1 ,x 2 ) = (0\Tij( Xl )ij(x 2 )\P). (4) 
where T represents time-order product, and transform it into the momentum space 

Xp(q) = e~ tP - x f d*xe-^ X (xi, x 2 ). (5) 



Here P is the four-momentum of the meson and q is the relative momentum of quark and 
antiquark. We use the standard center of mass and relative variables: 

X = 771X1 + 772X2, x = xx — x 2 , (6) 

where 77^ = m ^ m2 (i = 1,2). Then in momentum space the bound state BS equation reads 

{i>\ - m 1 ) X p{q)(fc + ma) = ^ J d*kG(P, q - k) X p(k), (7) 
where p\ and p 2 represent the momentum of quark and antiquark respectively. 

pi = r]iP + q, p 2 = r] 2 P - q. (8) 

G(P, q — k) is the interaction kernel which acts on x an d is determined by the interquark 
dynamics. Note in Eq.(|7]) mi and m 2 represent the effective constituent quark masses so that 
we could use the effective free propagators of quarks instead of the full propogators. This 
is an improtant approximation and simplification for light quarks. Furthermore, because 
of the lack of a fundamental description for the nonperturbative QCD dynamics, we have 
to make some approximations for the interaction kernel of quarks. In solving Eq.(0), we 
assume the kernel to be instantaneous (but with some retardation effect in a modified form 
for the kernel) and neglect the negative energy projectors in the quark propagators, because 



in general the negative energy projectors only contribute to quantities of higher orders due 
to M — E\ — E2 <C M + Ex + E2, where M, E±, and E2 are the meson mass, the quark kinetic 
energy, and the antiquark kinetic energy respectively. Based on the above assumptions 
the BS equation can be reduced to a three-dimensional equation, i.e., the reduced Salpeter 
equation, for the three dimensional BS wave function 

®p(q)= Jdq Xp(q°,q) } (9) 
(po _ Ei _ E 2 )§p(q) = A^ 7 ° J d 3 kG(P, q, k)§ p(k)^ A 2 _. (10) 

Here 

A^ = ^(£ 1 + 7 7 -Pi+™i7°), 

A 2 _ = J-(£ 2 - ^ . ^ _ m27 o ); (n) 

are the remaining positive energy projectors of the quark and antiquark respectively, and 
Ei = \Jm\ + pi 2 , E 2 = \Jrri 2 , + p~2 2 ' . The formal products of G$ in Eq.fllOD take the form 

G$ = Y l G i O&O i = G& + llt ®'fG v $, (12) 

i 

where = 7/J corresponding to the pertubative one-gluon-exchange interaction and = 1 
for the scalar confinement potential. 
From Eq.fllOD it is easy to see that 

A^p(q) = * p @, 

®p(q)A 2 _ = $p(q). (13) 

Considering the constraint of Eq . (|13"D , and the requirement of space reflection, in the rest 
frame of the meson (p = 0) the wave function <&p(q) for the 0~ and 1~ mesons can be 
written as 



$°r(g)=A^ 7 (i + 7 ) 757 A 2 ^(9), 

$L"(g) = Ah°(l + 7°) fr°^f(q), (14) 

where ^ = 7jU e M , e M is the polarization vector of 1~ meson, and (f(q), f(q) are scalar functions 
of (f . It is easy to show that Eq.(|I~4|) is the most genernal form for the 0~ and 1~ (S-wave) 



qiq 2 meson wave functions at the rest frame [e.g., for the 0~ meson wave function there are 
four independent scalar functions but with the constraint of Eq.(|i3|) those scalar functions 
can be reduced to one and expressed exactly as Eq.([L4])]. 

Substituting Eqs. (|T2|) and ( |Hj ) into Eq.(J10|), one derives the equations for </?(<?) and f{l) 
in the meson rest frame 1 



M Vx {q) = (E 1 +E 2 )< fl (q) 

E X E 2 + mim 2 + q 



where 
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Eim 2 + E 2 m\ 
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E x — E 2 
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r^s _ (mi + m 2 + Ey + E 2 )(E 1 m 2 + E 2 mi) 
AEiE 2 {m x + m 2 ) 



(E l +E 2 )f 1 (q) 

' r d 3 k(G s (q, k) - 2G v (q i k)){E im2 + E 2 m 1 )f x (k) 
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Mk), (15) 
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d 3 k(G s (q,k)+AG v (q,k))(q ■ k)fi(k) 
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fi(l) 
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fi(k), 



AEiE 2 



(17) 



(18) 



Eqs.flHf) and (|17D can also be formally expressed as: 



(M - E x - E 2 )tpi(q) = \d 3 k Ff r tf,1c)G i (?,iz)<t>i&), 

J i=S,V 

(M — Ei — E 2 )h{q) =Jd 3 kJ2 F f{q, k)Gi(q, k)fi(k). (19) 



i=S,V 



In most cases, the interaction kernel is of the convolution type, i.e., G(q,k) = G(q — k) = 
G(p), where p = q — k is the momentum exchanged between the quark and antiquark. In 
the nonrelativistic limit for both quark and antiquark, Eq. (plj) and fllTD can be expanded in 
terms of q /mi and q /m 2 , and they are identical with the Schrodinger equation to the 
zeroth order, and with the Breit equation to the first order. 

Ill Interaction kernel and retardation for confinement 

To solve Eq.(0) one must have a good command of the potential between two quarks. 
At present, the reliable information about the potential only comes from the lattice QCD 
result, which shows that the potential for a heavy quark-antiquark pair QQ in the static 
limit is well described by a long-ranged linear confining potential ( Lorentz scalar Vs ) and 
a short-ranged one gluon exchange potential ( Lorentz vector Vy ), i.e,|], |j, 

V s (r) = Ar, V v (r) = (20) 

6 r 

The lattice QCD result for the QQ potential is supported by the heavy quarkonium spec- 
troscopy including both spin-independent and spin- dependent effects|3|, f|, Here, as the 
first step, we will employ the static potential below regardless of whether the quarks are 
heavy or not 

V{r) = V s {r)+^®rVv{r\ 

V s (r) = Xr { - ?—L, 

ar 

Vv(r) = -l^M e — , (21) 
6 r 

where the introduction of the factor e~ ar is to avoid the infrared(IR) divergence and also to 
incorporate the color screening effects of the dynamical light quark pairs on the "quenched" 
QQ potential [|T2j| . It is clear that when ar<l the potentials given in (^1|) become identical 
with that given in (p0|). In momentum space the potentials are 



G(p) = Gs(P)+i^rG v (p), 
GsiP ' 



Gv(P] 



2 a s {p) 



(22) 



where a s (p) is the well known running coupling constant and is assumed to become a 
constant of 0(1) as V ^ 

12vr 1 

a s (p) = — — . (23) 

27 Ha + jf-) 

1V QCD 

The constants A, a, a, and Aqcd are the parameters that characterize the potential. 

Next, an important step is to take the retardation effect of scalar confinement into consid- 
eration. As disscussed in section I, the retardation effect of confinement will be approximatly 
treated by adding a retardation term ^ to the instantaneous part as given in Eq.@,and 
Pq will be treated to take its on-shell values which are obtained by assuming that the quarks 
are on their mass shells. Then this retardation term will become instantaneous (but not 
convoluted). This modified scalar confinement potential will include the retardation effect 
and become 



G s (p)^G s (p,k) = --5 3 (p)+ X ' 



a tt 2 (p 2 + a 2 ) 2 



2A 1 



+ —(^(p + k) 2 + m 2 - ^k 2 + m 2 ) 2 . (24) 

n 2 (p z + a 2 ) 6 

This shows that the retardation effect of confinement is taken into consideration in the 
way that the interaction kernel depends not only on p (the momentum exchanged between 
quark and antiquark) but also on k (the momentum of the quark itself). By calculation 
below one may see that the retardation effect is not negligible for cc states and become very 
significant for light-quark systems and we find it might be an useful remedy for the "intrinsic 
flaw" of the reduced BS equation with static scalar confinement. In the computation in the 
next section we will use 

A = 0.183GW 2 , a = 0.06GeV, a = e = 2.7183, A QCD = 0.15GeV, (25) 
All these numbers are in the scopes of customarily usage. 

IV Results and discussions 

Based on the formulism above, we have calculated the mass spectrum of quarkonium 
including both heavy- and light-quark systems. The numerical results with retardation are 
listed in Table I. The quark masses for the fit in Table I are 



m u = 0.35GeV, m d = 0.35GW, m s = 0.5GeV, 



m c = lMGeV, m b = A.83GeV. 



(26) 
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TABLE I. Calculated mass spectra of bb, cc, ss, and uu or dd states using reduced Salpeter 
equation with retardation for scalar confinement. The Experiment data are taken from Ref . fi~3|| . 
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TABLE II. The 2S — IS and 3S — IS energy level spacings of vector mesons with retardation 

and without retardation for scalar confinement 



For comparison with the results obtained without retardation, in Table II we give a list of 
25—15 and 35—15 energy level spacings for vector mesons in two cases, i.e., with retardation 
( using (f24"l)) and without retardation (using (p2|) ) for the scalar confinement potential. 



From Table I and Table II we can clearly see the following. 

(1) The calculated level spacings without retardation are generally smaller than their 
experimental values. This trend is already appreciable for charmonium and becomes a serious 



might improve the fit by readjusting the parameters (e.g., by enlarging the string tension), 
and this may work for low-lying heavy quarkonium states (e.g., cc states) but can not give 
a good global fit for high-lying states especially for light quarkonium states. 

(2) By adding the retardation term to scalar confinement potential the calculated level 
spacings are significantly improved. The fit for bb and cc states is very good, and the fit for 
light vector mesons is also good, while the fit for light pseudoscalar mesons such as the pion 
is poor, which is probably due to the fact that the light pseudoscalar mesons are essentially 
Goldstone bosons and therefore the instantaneous and on-shell approximations no longer 
work well for them. 

As emphasized in Section I, the approximate treatment for the retardation effect (in par- 
ticular, the on-shell approximation) of scalar confinement should be good for heavy quarko- 
nium states. Indeed, it has been shown Q that for heavy quarkonium in nonrelativistic ex- 
pansion the role of the retardation is just to cancel the troublesome term, -^(j^S 1 + Sp 2 ) 
in Eq.(^) and then remove the bad disturbance to the mass spectrum. 

For light quarkonium states with constituent quark masses m u = ~ SSOMeV^the 
retardation effect bocomes even more significant. In these systems nonrelativistic expansion 
is no longer good, but we can see the physical effect of retardation through an extreme case. 
That is the zero quark mass limit, which has been used for analysing the "intrinsic flaw" of 
scalar confinement, indicating that light quarks can only have very weak confinement if it 
is an instantaneous Lorentz scalar potential ||. To see how the "retardation" part changes 
the trend of light quarks seeing a weaker confining potential than heavy quarks at large 
distances, it is useful to consider again the limit of zero quark mass. As light quark systems 
are more sensitive than heavy quark systems to the behavior of interaction at large distance, 
we will restrict our discussion only on the scalar confinement potential part. 

In the zero-quark mass limit, as m q — > 0, the coefficients for the scalar potential G s in 
Eq.(|l9D for the 0~ and 1~ mesons will reduce to 




problem for light quarkonium states. This result agrees with that obtained in Ref. 0. We 





1 



q ■ k 



— * — * 

where q = \q\, k = \k\. It is clear that these coefficients will vanish when q — > k. On 



the other hand, however, the static linear confining potential in momentum space behaves 
as G s (q — k) oc (q — k)~ 4 and is strongly weighted as q — > k in Eq. (^) . This is the 
reason why the light quarks can only have weak confinement, which leads to very narrow 
energy level spacings for light quarkonium states. This bad situation will be changed if the 
retardation is taken into account. In fact, the covariant form of confinement interaction may 
take the form Gs{q,k) oc [(q — k) 2 — (q — k ) 2 }~ 2 and in the on-shell approximation that 
go = m 2 + (f-, k\ = m 2 + k 2 it becomes 

G s (q, k) oc [-2m 2 + m 2 - + m 2 \ + 2(qk - q- k)]- 2 , (28) 

q k 

where in the zero-quark mass limit p, k ^> m — > 0. We can see immediately that it has 
two distinct features from the static confining potential Gs{q, k) oc [(q — k) 2 }~ 2 . First, with 
retardation the scalar interaction Gs{q, k) is strongly weighted when q and k are co-linear 
(q || k) that leads to q- k — qk, whereas the static linear potential only peaks at q = k. This 
indicates that the former is strongly weighted in a much wider kinematic region than the 
latter. Second, when q || k, Gs{q, k) oc 0(m~ 4 ), this mass dependence, which is absent in the 
static linear potential, will enhance confinement interaction and overwhelm the suppression 
factor apppearing in the coeficients F$ and Fg . As a result, even in the zero-quark- 
mass limit the effective scalar interaction will not be weakened, and this is just due to the 
retardation effect in the scalar interaction. 

In practice, for the constituent quark model, which is essentially used in the present work, 
the quark mass can not be zero, and the on-shell approximation may not be as good as for 
heavy quark systems. But in any case the analysis given above in the zero mass limit is useful 
for understanding the qualitative feature of the retardation effect in quark confinement. 

In this paper, we have tried to clarify the problem pointed out by Durand et al. for the 
static scalar confinement in reduced Salpeter equation. The " intrinsic flaw" of the Salpeter 
equation with static scalar confinement could be remedied to some extent by taking the 
retardation effect of the confinement into consideration. In the on-shell approximation for 
the retardation term of linear confinement, the notorious trend of narrow level spacings for 
quarkonium states especially for light quarkonium states is found to be removed. A good fit 
for mass spectrum of S-wave heavy and light quarkonium states (except the light pseudoscalar 
mesons) is obtained using one-gluon exchange potential and the scalar linear confinement 
potential with retardation taken into account. Although for light quark systems the on-shell 
appproximation may not be good, the qualitative feature of the retardation effect is still 
manifest. We may then conclude that at phenomenological level including the retardation 
effect into the scalar confinement may be necessary and significant. Nevertheless, it is still 



premature to assess whether or not quark confinement is really represented by the scalar 
exchange of the form of (p 2 — pl)~ 2 , as suggested by some authors as the dressed gluon 
propogator to implement quark confinement. We hope that our investigation can provide 
some useful information for the understanding of confinement. Further discussions concening 
heavy-light mesons will be given in another publication. 
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